arXiv:1508.04861vl [nucl-th] 20 Aug 2015 


EPJ manuscript No. 

(will be inserted by the editor) 


Hyperon Puzzle, Hadron-Quark Crossover 
and Massive Neutron Stars 

Kota Masuda^’^ Tetsuo Hatsuda^’^ and Tatsuyuki Takatsuka^ 

^ Department of Physics, The University of Tokyo, Tokyo 113-0033, Japan 
^ Theoretical Research Division, Nishina Center, RIKEN, Wako 351-0198, Japan 
^ Kavli IPMU (WPI), The University of Tokyo, Chiba 277-8583, Japan 

Received: date / Revised version: date 

Abstract. Bulk properties of cold and hot neutron stars are studied on the basis of the hadron-quark 
crossover picture where a smooth transition from the hadronic phase to the quark phase takes place at 
finite baryon density. By using a phenomenological equation of state (EOS) “CRover” which interpolates 
the two phases at around 3 times the nuclear matter density (po), it is found that the cold NSs with the 
gravitational mass larger than 2Mq can be sustained. This is in sharp contrast to the case of the first-order 
hadron-quark transition. The radii of the cold NSs with the CRover EOS are in the narrow range (12.5T0.5) 
km which is insensitive to the NS masses. Due to the stiffening of the EOS induced by the hadron-quark 
crossover, the central density of the NSs is at most 4 po and the hyperon-mixing barely occurs inside the NS 
core. This constitutes a solution of the long-standing hyperon puzzle. The effect of color superconductivity 
(CSC) on the NS structures is also examined with the hadron-quark crossover. Eor the typical strength 
of the diquark attraction, a slight softening of the EOS due to two-flavor CSC (2SC) takes place and the 
maximum mass is reduced by about 0.2 Mq. The CRover EOS is generalized to the supernova matter at 
finite temperature to describe the hot NSs at birth. The hadron-quark crossover is found to decrease the 
central temperature of the hot NSs under isentropic condition. The gravitational energy release and the 
spin-up rate during the contraction from the hot NS to the cold NS are also estimated. 

PACS. 21.65.Qr Quark matter - 26.60.-c Nuclear matter aspects of neutron stars - 97.60.Jd Neutron 
Stars 


1 Introduction 

Strongly interacting matter described by quantum chro¬ 
modynamics (QCD) is believed to have a rich phase struc¬ 
ture under the change of external parameters such as the 
temperature (T) and the baryon chemical potential (/r) 
[T]. At low T and /i, the system is in the hadronic phase 
where the dynamical breaking of chiral symmetry and con¬ 
finement of quarks and gluons take place. At high T and 
low the quark-gluon plasma with deconfined quarks and 
gluons was predicted theoretically and is under active ex¬ 
perimental studies by the relativistic heavy-ion collisions 
at RHIC and LHC [5]. At low T and high fi, the supercon¬ 
ducting quark matter with deconfined quarks is expected 
to appear, which is relevant to the central core of neutron 
stars [3]. 

The transition from the hadronic matter to the quark- 
gluon plasma at high T has been studied quantitatively 
by using the lattice QCD simulations [J]. On the other 
hand, the transition from the hadronic matter to the quark 
matter at high p, is not well understood partly due to 
the lack of reliable first-principle theoretical methods; the 
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Monte Carlo simulations in lattice QCD are not suitable 
for fi/T ^ 1 because of the fermion sign problem [T]. 

Under such circumstance, any information from neu¬ 
tron stars (NSs), whose central cores may reach the baryon 
density relevant to the hadron-quark transition, is quite 
useful to understand the structure of high density matter. 
Among various observables for NSs [5], the surface temper¬ 
ature (Ts) the mass (M), the radius (R) and the magnetic 
field {B) are particularly informative. Indeed, a massive 
NS (PSR J1614-2230) with M = (1.97±0.04)Mo observed 
through the Shapiro delay technique [S] and another mas¬ 
sive NS (PSR J0348-b0432) with M = (2.01 ± 0.04)Mo 
[7] may give stringent constraints on the equation of state 
(EOS) of dense matter. 

Historically, Gibbs phase equilibrium conditions have 
been adopted to describe the transition between the hadronic 
matter composed of point-like hadrons and the quark mat¬ 
ter composed of weakly interacting quarks. However, in 
the transition region, neither the assumption of point-like 
hadrons nor that of weakly-interacting quarks are justi¬ 
fied. Indeed, there may arise a smooth crossover between 
the hadronic matter and the quark matter: Furthermore, 
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Fig. 1. Schematic pictures of the QCD phases as a function of the baryon density (p) under the assumption of (a) the first-order 
transition and (b) the hadron-quark crossover. The mixed-phase region in (a) and the crossover region in (b) are indicated by 
the shaded area. 


the system in the crossover region would be strongly in¬ 
teracting [ 8 ]. 

Fig. 0 illustrates the difference between (a) the first- 
order transition where the phase separation between 
hadrons and quarks takes place, and (b) the crossover 
where the percolation of finite size hadrons takes place. 
Such a percolation picture of hadrons has been first dis¬ 
cussed in Refs. wm. and later elaborated in the contexts 
of the hadron-quark continuity [ma and the hadron- 
quark crossover UMM- 

Recently, the present authors have shown that the 
hadron-quark crossover at around three times the normal 
nuclear matter density (po = 0.17 fm“^) can lead to a 
stiffening of EOS and sustain the 2Mq NSs [T5llT6lfT7] in 
contrast to the case of hadron-quark first-order transition. 
Also, it was shown that such a stiffening due to hadron- 
quark crossover can avoid the so-called “Hyperon Puzzle” 
as discussed in §2. (See also the related works [rSlITOll^ .i 

In this article, we discuss bulk properties of cold and 
hot NSs on the basis of the new EOS with the hadron- 
quark crossover (the “CRover” EOS) introduced in [ISlITHl 
[H]. In §2, we summarize the conventional hadronic EOS 
with and without hyperons. In §3, we summarize detailed 
properties of hadronic EOS to be used throughout the 
present study. In §4, we summarize the quark EOS based 
on the (2-|-l)-flavor Nambu-Jona-Lasinio (NJL) model at 
high density. In §5, We introduce a phenomenological ap¬ 
proach to interpolate the hadronic matter and the quark 
matter. In § 6 , we show the bulk properties of cold neutron 
stars using the CRover EOS at T = 0 (abbreviated as cold 
CRover EOS) with and without color superconductivity. 
In §7, we show the bulk properties of hot neutron stars 
at birth using CRover EOS at T 7 ^ 0 (abbreviated as hot 
CRover EOS). §8 is devoted to summary and concluding 
remarks. 


2 Hyperon Puzzle 

Let us first consider what would be the neutron star struc¬ 
ture under the hadronic EOS with and without hyperons 
(Y=A, IJ~, S^). Although there exist large uncertain¬ 
ties for the two-body Y-N interactionsvarious differ¬ 
ent models suggest that hyperons may appear at densities 
of several times po- Those hyperons introduce signihcant 
softening of EOS and even the well-established 1.4Mq NSs 
may not be reproduced (see the reviews, mm and the 
references therein.) The three-body force in the hyperon 
sector originally suggested in may or may not describe 
the 2Mq NSs depending on its strength pniMl[^[l26l[?7] . 
This is called the “Hyperon Puzzle”. 

To see the “Hyperon Puzzle” more explicitly, let us 
take four hadronic EOS with hyperons, TNI2u, TNI2, 
AVIS-bTHE-bAA, and SCLSAr. Here TNI2u (TNI2) is 
the EOS based on the non-relativistic G-matrix approach 
with the incompressibility k = 250MeV and with (with¬ 
out) the hyperon three-body force. AV18-|-TBF-|-ylA' [32] 
is also based on the non-relativistic G-matrix approach 
with the AVI 8 nucleon-nucleon potential, the Urbana- 
type three-body nucleon potential and the Nijmegen soft¬ 
core nucleon-hyperon potential. SCL3AA' [53| is a rela¬ 
tivistic mean-held model with chiral SU(3) symmetry. As 
a typical nuclear EOS without the hyperons, we take APR 

O' 

In Eigj^a), we plot the hadronic EOS with hyperons 
(TNI2u, TNI2, AV18-kTBE-kAA and SCL3Ar) together 
with APR. Pilled circles on each line denote the density 
where the hyperon-mixing starts to occur. One can see 
that (i) the mixture of hyperons softens the equation of 
state relative to APR, and (ii) onset of the hyperon-mixing 
is shifted to higher density if we consider the three-body 
interaction in the hyperon sector. In Figj^b), the M — R 

^ We note that it is important to have new data on hyperon 
interactions from the Y-N scattering and hypernuclei exper¬ 
iments at J-PARC |28ll29il^ and also from the lattice QCD 
simulations at the physical quark masses m- 
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Fig. 2. (a) The hadronic equation of states with and without the hyperon mixing. The black line denotes the EOS without 
hyperon, APR (AV18-|-5w-|- UIX*). The red lines correspond to the EOS with hyperons; TNI2 (only three-nucleon force with 
K = 250MeV) and TNI2u (universal three-baryon force with k = 250MeV). The blue line corresponds to AV18-|-TBF-|-ylX' 
(G-matrix with AV18 -I- 3-nucleon force -I- hyperons), and the green line corresponds to SCLSAX' (relativistic mean-field model 
with chiral SU(3) symmetry), (b) M — R relations for five EOSs considered in (a). In both figures, filled circles on each line 
show the density where the hyperons start to mix. Also, the gray band denotes M = (1.97 ± O.O4)M0 for PSR J1614-2230 [B] 
and the solid horizontal line shows M = 1.44M0 corresponding to PSR 1913-1-16. Those figures are adapted and modified from 

mM- 


relations with these EOSs are plotted by the same color 
lines. The gray band shows M = (1.97 ± O.O4)M0 corre¬ 
sponding to PSR J1614-2230 [5] and the solid horizontal 
line shows M = 1.44M0 corresponding to PSR 1913-1-16. 
In Table [l] we summarize the nuclear incompressibility 
K, the threshold density for A and E~, the maximum 
mass Mmax, the radius and the central density pc for these 
hadronic EOSs with hyperons. EOSs become soft drasti¬ 
cally due to the emergence of hyperons. 

3 Hadronic EOS with Hyperons 

Since we will take TN12u and TN12 in the following anal¬ 
yses, we summarize here how these hadronic EOS with 
hyperons are obtained 

1 . Effective two-baryon potentials Vbb' {B = n, p, A, 
S~ ) are constructed on the basis of the G-matrix for¬ 
malism to take into account their density-dependence. 

2. A phenomenological thee-nucleon interaction expressed 
in a form of two-body potential Unn' |3S] is introduced 
to reproduce the saturation of symmetric nuclear mat¬ 
ter (the saturation density po = 0.17 fm“^ and the 

2 “TNI” implies that the Three-Nucleon Interaction is taken 
into account, “2” implies k = 250 MeV, and “u” implies 
that the three-body interaction is introduced universally in the 
octet baryon sector. 


binding energy Eq = —16 MeV) and the incompress¬ 
ibility K = 250 MeV compatible with experiments. 

3. Universal three-body repulsion is assumed for the hy¬ 
perons in TNI2u through the replacement, Unn' —t 
Ubb', which is necessary to sustain 1.4M0 even with 

Table 1. Properties of various hadronic EOSs with hyperons; 
TNI2, TNI2u, AV18-|-TBF-|-ylX' and SCL3/IV. k is the nuclear 
incompressibility and pth is the threshold density of hyperon- 
mixing with po being the normal nuclear density. R and pc 
denote the radius and central density of the maximum mass 
(Muiax) NS, respectively. The numbers in the parentheses are 
those without hyperons. “*”s indicate that the numbers are 
read from the figures in m- 


EOS 

TNI2 

TNI2u 

AV18-bTBE 

SCL3Ar 

K (MeV) 

250 

250 

192 

211 

Pth(/l)/po 

2.95 

4.01 

2.8* 

2.24 

pth(V")/po 

2.83 

4.06 

1.8* 

2.24 

Mmax /Mq 

1.08 

1.52 

1.22 

1.36 


(1.62) 


(2.00) 

(1.65) 

R(km) 

7.70 

8.43 

10.46 

11.42 


(8.64) 


(10.54) 

(10.79) 

Pc/ Po 

16.10 

11.06 

7.35 

6.09 

(9.97) 


(6.53) 

(6.85) 
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(a) (b) 

Fig. 3. (a) The number fractions {nu,d,s,e/ntot with ntot = Uu + nd + ris = 3p) as a function of the baryon density p. Red, blue 
and green lines correspond to u, d and s quark, respectively. The black line corresponds to the electron number fraction x 100. 
Muon does not appear due to the emergence of s quarks, (b) The constituent quark masses {Mi) as a function of p. Colors on 
each line are the same with those in (a). These figures are adapted from EH. 


hyperons, while the three-body repulsion is introduced 
only in the nucleon sector in TNI2. 

4. By using Vbb' + Ubb', we calculate the particle com¬ 
position yi {i = n, p, A, S~, e~ and pT) under charge 
neutrality and /3-equilibrium to obtain the hadronic 
EOS as a function of total baryon density p at T = 0. 

As shown in Fig|^a), TNI2u EOS is moderately stiff even 
with hyperon-mixing, but the corresponding maximum 
mass of NS is I.SMq, so that it is not sufficient to re¬ 
produce 2Mq NSs as shown in Figj^b). 

4 (2+l)-flavor Quark EOS with color 
superconductivity 

At high baryon density relevant to the central core of the 
NSs, baryons are started to overlap with each other and 
quark degrees of freedom may arise with p ~ 400 — 500 
MeV. However, at such chemical potentials, the QCD run¬ 
ning coupling is still large and the quark matter would be 
strongly interacting. Analogous strongly interacting mat¬ 
ter at finite T was originally discussed theoretically in [SZl 
[381IM] and was later studied experimentally in the rela¬ 
tivistic heavy-ion collisions at RHIC and LHC P]- 

Under the situation that the Monte Carlo simulations 
in lattice QCD are not applicable to p/T ^ 1 due to the 
sign problem, we adopt the (2-1-1)-flavor Nambu—Jona- 
Lasinio (NJL) model which captures the essential prop¬ 
erties of QCD such as the dynamical breaking of chiral 
symmetry and its restoration at finite T and p goEl]. 


The model Lagrangian we consider reads 
1 ® 

>Cnjl = qii0 - m)q + -Gg ^[{qX°-q)^ + {qi"t5X°'qf] 

a—0 

-Go [detg(l -b 75)9 + h.c.] - -g^ {qi^qf 
+ y X! (9*75A^r^Gg'^)(g^Gf75A^T^9) 

^ /,A^2,5,7 

[{qKTk{i + ib)cf) 

{qXjTu{l + 75)Gq)(gi(l -b 75)gj) -b h.c.] 

where the quark field qf has three colors labeled by a and 
three flavors labeled by i with the current quark masses 
mi. The term proportional to Gg is a t/(3)z, x U{3)ii sym¬ 
metric four-fermi interaction where A“ are the Gell-Mann 
matrices in flavor space with A® = \/2j3 I. The term pro¬ 
portional to Gg, is called as the Kobayashi—Maskawa—’t 
Hooft (KMT) coupling which breaks U{1)a symmetry ex¬ 
plicitly The term proportional to g^{> 0) gives a 

universal repulsive force among different flavors. The term 
proportional to H gives a diquark condensation with color 
anti-triplet, flavor anti-triplet and = O’*" channel where 
G = 17270 is the charge conjugation matrix and r“ are 
the Gell-Mann matrices in color space with = y^2/3 I. 
The term proportional to G'^ is obtained by the Fierz 
transformation of the KMT term and gives a coupling be¬ 
tween the chiral condensate and the diquark condensate. 
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(a) (b) 

Fig. 4. (a) The gap parameters Ai {i = 1, 2,3) as a function of the quark chemical potential /i. Solid red line: A 3 corresponding to 
the pairing between u quark and d quark. Solid blue line: Zii and A 2 corresponding to the ds pairing and su pairing, respectively, 
(b) The energy density (e) and pressure (P) as a function of baryon number density p. The red lines show the quark EOSs with 
diquark condensate. The blue lines show the quark EOSs without diquark condensate. 


Here we have introduced a notation, Ai^ 2.3 = A 7 _ 5_2 and 

^ 1 , 2,3 = '^7,5,2- 

In the mean field approximation, the constituent quark 
masses Mi and the gap parameters Ai are generated dy¬ 
namically through the NJL interactions. 


G' 

Mi = rrii — 2Gg(Ti -|- 2G^ajak H—(1) 


A = - 



Si, 


( 2 ) 


where ai = (qiqi) is the chiral condensate in each flavor, 
Si = {q^C^^Xifiq) is the diquark condensate in each color 
and flavor with corresponding to the cyclic per¬ 

mutation of u, d and s. The thermodynamic potential 12 
is related to the pressure as 12 = —TlogZ = —PV with P 
given by 


P{T,pLu,d,s) 


T 

I 


E 


Trlnl 

A5 i(wf,p)\ 

/ (27r)3 ' 

V T ) 


G g ^ ^ CT I 4G 
i 

- E liH-G'J\s,f. 

z=l,2,3 



is the quark propagator in the Nambu-Gor’kov represen¬ 
tation, 

fc-ilb _ / [^ 0 "] ^ Si=i,2.3 \ 

where 

[G±]-i=7^-M±7oA‘=®. (4) 

Here, M is a unit matrix in color space and a diagonal 
matrix in flavor space, diag(M„, M^, Mg). The effective 
chemical potential matrix is defined from 

A®® = A - 5v E '^3 

3 

where each component of jl is given by 

/^a/3 = + ^^QQ"^)Sa0 + (/i3(T-3)a/3 + A<8('r8)a/3)d*G 

There are nine independent parameters in the (2-1-1)- 
flavor NJL model; the UV cutoff, A, the coupling con¬ 
stants, Gg, G^ ,gy,H and G'^ and the current quark masses, 
mu,d and mg. Five parameters except for gy,H and G'^ 
have been determined from hadron phenomenology in the 
vacuum. In this article, we adopt the HK parameter set 
[4fi] (results for other parameter sets, see [IS]); 


Here ii/e = {2£ + 1 )ttT is the Matsubara frequency, rii = A = 631.4MeV, G^A^ = 3.67, G^A^ = 9.29, 

{qlqi) is the quark number density in each flavor, and S rriu^d = 5.5MeV, mg = 135.7MeV. (6) 
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(a) 



p/p^ 


Cb) 


Fig. 5. (a) The relation between the interpolated energy density and the resultant pressure for = 0 (green) and = 0.5G^ 
(blue) without CSC, H = G'jy — 0. The crossover window is (2 — 4)po- The filled circle denotes the threshold density of 
the strangeness, (b) Sound velocity squared Vs as a function of baryon density p. Solid lines show Vs — p obtained from the 
interpolated EOS with Qy jG^ = 0, 0.5, while the dotted line corresponds to the TNI2u hadronic EOS. The filled circles denote 
the points beyond which strangeness starts to appear. Those figures are adapted and modified from m- 


For Qy , we change its magnitude in the following range 

mi^, 

0 < If < 0.5, (7) 

The parameters H and G'^ are chosen to be 
and G’jj = G^ as characteristic values. (Corresponding 
values from the Fierz transformation are H = and 
G'jj = Gj 3 . For extensive analyses with other choice of 
parameters in the diquark channels, see [46llT7ll48] .) 

The EOS of quark matter with strangeness is obtained 
from the above model under two conditions: (i) the charge 
neutrality among quarks and leptons, ‘^riu — ^rid — |ns — 
ne — = 0, (ii) the color neutrality among quarks, = 

ng = rif,, and (iii) the /3-equilibrium among quarks and 
leptons, Pd = Ps = Pu + Pe and pe = Pg,- The ground 
state of the system is obtained by searching the maximum 
of the pressure in Eq. ([^ with the conditions, 

dP dP dP 
7)-= TFa-= 7)-= 

OCru,d,s 0^1, 2 ,3 Op3^s 

Let US first discuss a composition of the /3-equilibrated 
quark matter at T = 0 without color superconductiv¬ 
ity {H = G'^ = 0). In FigjS (a), the number fractions 
{nu,d,s,e/ntot with ritot = Uu + nd + Us = 3p) as a func¬ 
tion of the baryon density p are plotted. In Eigj^(b), the 
constituent quark masses (Af) as a function of p are plot¬ 
ted. These figures do not depend on the magnitude of the 
vector type interaction py because the flavor-independent 


( 7 ^,-type interaction leads to a pressure in Eq. § depend¬ 
ing only on 

At low densities, the s quark appears only above pth — 
4po due to its heaviness as can be seen from Figj^(a): Here 
Pth is determined by the condition, ps{pth) = Afs(pth)- 
The dynamical masses of u and d quarks decrease rapidly 
below Pth due to partial restoration of chiral symmetry, 
while the s quark is affected only a little through the KMT 
interaction as seen from Figj^(b). Once the s-quark whose 
electric charge is negative starts to appear above pth, the 
number of electrons decreases to satisfy the charge neu¬ 
trality. Since the electron chemical potential does not ex¬ 
ceed the muon mass 106MeV, the muons do not appear 
even at high density. In the high density limit, the system 
approaches to the flavor symmetric u, d, s matter without 
leptons. 

Once the interactions in the diquark channels {H and 
G'^) are switched on, the color superconductivity (CSC) 
develops as shown in Fig. |^(a) where the gap parameters 
Ai as a function of the quark chemical potential p with 
9vl^s = 0.5 are plotted. The red line corresponds to the 
ud pairing, and blue line corresponds to the ds or su pair¬ 
ings. With H/Gg = 1, two-color superconductivity (2SC) 
appears as soon as the baryon density of the quark mat¬ 
ter becomes finite at p = 335 MeV. Then the first-order 
transition from the 2SC to the color-flavor-locking (CFL) 
phase takes place at around p = 520MeV. The diquark 
condensates affect the behavior of the number fractions 
{rii^e/ntot and the constituent quark masses {Mi) through 
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(a) (b) 


Fig. 6. (a,b) Comparison between the CRover EOS and the sound velocity squared without CSC {H = = 0) and those 

with CSC {H = Gs and G'd = Go) for g^/G^ = 0.5. 


the coupled equations, Eq.(|^: Those are taken into ac¬ 
count into the following calculations with CSC. 

Fig. 13(b) shows the thermodynamic quantities {P and 
e) as a function of p in the (2-|-l)-flavor NJL model. The 
red (blue) lines correspond to the case with (without) 
CSC. In terms of the baryon density, 2SC (CFL) appears 
for 2po < p < 5po (p > 5po) in the present parameter set. 
The plateau of the red line {P with CSC) reflects the fact 
that there is a first-order transition from 2SC to CFL. As 
we will see later, baryon density relevant to neutron stars 
with the hadron-quark crossover is below 5po- Therefore, 
the CFL phase barely appears in the central core of the 
NSs in the present model. 


5 CRover: New EOS with Hadron-Quark 
Crossover 

We now introduce the following phenomenological inter¬ 
polation of the energy per baryon E = E/N^ at T = 0 

pfolfTH] 


E{p) = Eb_{p)w-{p) + Eq{p)w+{p), (9) 

where Fh and Eq represent energy per baryon in the 
hadronic matter with leptons and that in the quark mat¬ 
ter with leptons, respectively. W- and = 1 — W- are 
the weight functions 


w±{p) = 


1 

2 


^litanh (^^4^)) ’ 


( 10 ) 


where p and P are the phenomenological parameters which 
characterize the averaged crossover density and the width 
of the crossover window, respectively. Similar weight func¬ 
tion has been previously used to interpolate the hadronic 
phase and the quark-gluon plasma at finite temperature 
snuns]. The window p — P^p^p + P characterizes the 
crossover region in which both hadrons and quarks are 
strongly interacting, so that neither pure hadronic FOS 
nor pure quark FOS are reliable. 

The other observables can be obtained by using the 
thermodynamic relations from E(p) = e/p; 


P = p^ 


vl{p) = 


,dE 
dp ’ 
dP 


de 


P = 


de 


dp' 
K 

e + P' 


K = p 


dP 

w 


( 11 ) 


Here K is the bulk modulus which must be positive semi- 
definite for thermodynamic stability. Also, Vg is the sound 
velocity which is a characteristic measure of the stiffness 
of the FOS. 

In Fig[^ (a), we show the interpolated FOS at T = 0 
(cold CRover FOS) with the TNI2u as a hadronic FOS 
and the NJL with = 0 [g^ = O.SGg) as a quark FOS. 
The diquark condensates are switched off here (iL = = 

0). The sound velocity squared as a a function of p is shown 
in Figj^(b). In both figures, the onset of the strangeness is 
indicated by the filled circles. From both figures, one finds 
that the CRover FOS becomes stiffer than the quark phase 
and the hadronic phase in the crossover region indicated 
by the shaded band. Such stiffening is induced by the ex¬ 
tra pressure originating from the derivative of p acting on 
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Fig. 7. (a) The M — pc relations with the CRover EOS (solid lines) and the TNI2u hadronic EOS (the dotted line). The 

crossover window is fixed to be (p, T) = (3po,Po)- The cross symbols denote the points of Mmax, while the filled circles denote 
the points beyond which the strangeness appears. The gray band denotes M = (1.97 ± O.O4)M0 for PSR J1614-2230. The solid 
black line denotes M — 1.44Mq for PSR 19134-16. (b) The M — R relation with the same EOSs as (a). These figures are adapted 
from [16] . 


w±{p) in Eq. 

P{p) = Pii{p)w-{p) + Pq{p)w+{p) + AP(p). (12) 

The extra pressure AP is required from thermodynamic 
consistency and has a property, AP{p = 0) = AP{p = 
oo) = 0 by definition, i.e. it is a function localized in the 
crossover region. 

By turning on CSC with p[ = Gs and G'jj = Go, the 
interpolated EOS becomes a little bit softer than the case 
without CSC in the crossover region as shown by the red 
lines in Fig|^(a) and (b). Associated with this, the onset 
density of the strangeness is reduced from 4p to 3.6po- As 
we have discussed in §4, there is little room for the CEL 
phase to appear inside NSs in our CRover EOS, it is not 
considered in this figure. 

6 Neutron Stars with CRover EOS at T = 0 

We study the structure of the spherically symmetric neu¬ 
tron stars in hydrostatic equilibrium by solving the Tolman- 
Oppenheimer-Volkov (TOV) equation; 

dP G 

— = -- (M(r) -f AirPr^) (e + P) {I - 2GM{r)/r)~^ , 

dr 

M{r) = f 47rr'^£(r')dr'. (13) 

Jo 

Here r being the radial distance from the center and G is 
the gravitational constant. 


6.1 Case without Color Superconductivity 

First, we consider the case without the diquark condensate 
(id = G'^ = 0) as a baseline. In Fig[^ (a), we plot M — 
Pc relation from the CRover EOS (interpolation between 
the TNI2u hadronic EOS and the NJL quark EOS with 
9vl^s — OjO.5 in the crossover region {p,r) = (3po,Po))- 
For comparison, the M — pc relation only with TNI2u 
hadronic EOS is plotted by the dotted line. Fig[^(b) shows 
the corresponding M — R relation for the same EOSs as 
Fig[^(a). Strong correlation between the sudden stiffness 
of the sound velocity in Figj^ (b) and the NS masses in 
Fig0 (a,b) can be seen. Also, as increases from zero, 
the quark EOS and hence the CRover EOS become stiffer, 
which increases the maximum mass (indicated by the cross 
symbol) beyond 2Mq. 

One important aspect of the present result is that the 
radii of NSs are confined in a narrow range: For example, 
in the case of Py/G^ = 0(0.5), all the NSs with 0.5 < 
M/Mq < 2.0 (0.5 < M/Mq < 2.5) have the radius in 
the range R = (12.0 ± 0.2) km (i? = (12.5 ± 0.5) km). 
Such a narrow window of R independent of the values of 
M will confront the present and future observations of the 
neutron star radii 

Another important aspects of the present result is the 
onset of the strangeness inside the NS core as indicated 
by the filled circles in Figj^b): The strangeness appears 
only in very massive NSs with nearly the maximum mass 
if we have hadron-quark crossover. In fact, pc = 2.4po for 
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(a) (b) 


Fig. 8. (a) The M — relations with the CRover EOS (solid lines) and TNI2u hadronic EOS (dotted line) with and without 
the CSC phase for /Gs = 0.5. (b) The M — R relation with the same EOSs as (a). 


2Mq while pc = 4.4/9o for 2.59Mq under the CRover EOS 
with = O.SGg, so that even the observed 2Mq NSs are 
unlikely to have strangeness inside. 


6.2 Case with Color Superconductivity 

In Fig. 1^ (a,b), M — pc and M — R relations are plot¬ 
ted by using the CRover EOS with and without CSC 
given in Figjm a). For comparison, the results of the TNI2u 
hadronic EOS are shown by the black dotted lines. As we 
discussed in the last section, the CSC softens the EOS. 
Then, the Mmax of the NS with CSC becomes smaller by 
0.2Mq than that without CSC phase. Such a reduction 
of Mijiax due to CSC is consistent with other calculations 
(see e.g., [S31[51])- 

Two remarks are in order here about the effect of CSC 
on the M — R relation: (i) The central density of the NSs 
does not exceed 4.9po in CRover EOS with CSC, so that 
the CEL phase barely appears inside the star, (ii) The 
effect of 2SC already becomes visible for low mass stars 
(M < 0.5Mq) with the central density below 2 pQ. This 
is because we have a smooth interpolation between the 
hadronic EOS and quark EOS, so that the 2SC compo¬ 


nent has small but non-negligible contribution even below 
2 pq. Physically, this could be interpreted as partial per¬ 
colation of the nucleons into quarks with strong diquark 
correlations. 

Table is a summary of the the maximum mass M^ax 
and the central density pc normalized by po with the 
CRover EOS. Two sets of hadronic EOS (TNI2u and TNI2) 
are adopted, but the difference is small. The strength of 
the repulsive vector interaction and the crossover den¬ 
sity p are changed to see the sensitivity of the results. As 
g^ becomes larger and the p becomes smaller, the max¬ 
imum mass Mmax increases due to the presence of the 
quark matter. The effect of CSC generally decreases Mmax 
and increases Pc- As long as p is around 3poj the CRover 
EOS can easily accommodate the 2Mq NSs. 

7 Neutron Star with CRover EOS at T 7 ^ 0 

In this section, we describe the properties of hot neutron 
stars created after the core-collapsed Type-II supernova 
explosion by generalizing the idea of the hadron-quark 
crossover discussed in the previous sections (For detailed, 
see [17]). 


Table 2. M^bx/Mq (pc/po) for different choices of hadronic EOS, quark EOS and the crossover window. 



flv = 0 

without CSC 

9v = 

without CSC 

Qv = U.5Gs 

with CSC 

(p,c) 

(3po, Po) 

(5po, 2po) 

(3po,po) 

(5po, 2po) 

(3po,po) 

TNI2u 

TNI2 

2.02 (4.5) 
2.02 (5.8) 

1.86 (8.7) 
1.84 (9.1) 

2.59 (4.4) 
2.59 (4.3) 

2.25 (6.1) 
2.23 (6.8) 

2.40 (4.9) 
2.40 (4.8) 


























10 


Masuda, Hatsuda and Takatsuka: Hyperon Puzzle, Hadron-Quark Crossover and Massive Neutron Stars 


BPS HS TNI2 


NJL 


10-3 1 2 2.4 p/po 



BPS HS TNI2_^ 



Fig. 9. Comparison between the hot NS and cold NS with M = 2Mq obtained by the CRover EOS with = O.SGs. We take 
{Yi,S) = (0.3,1) to characterize the hot NS. As for the details of the EOSs adopted at different densities, BPS, HS, TNI2 and 
NJL, see the text. This figure is adapted from m- 


During the first few seconds after the core bounce, the 
proto-neutron star with the radius ~ 200 km is formed. 
It undergoes a rapid contraction and evolves into a “hot” 
neutron star with the radius ^ 20 km in about 1 second (or 
to a black hole). The hot NS at birth in quasi-hydrostatic 
equilibrium is characterized by the supernova matter with 
the lepton fraction, Yi = Yg + Y,y ^ 0.3 — 0.4, and the 
entropy per baryon, S/N = S' ~ 1 — 2. The neutrinos 
are trapped inside the hot NSs and contribute to the /3- 
equilibrium. With this as an initial condition, the hot NS 
contracts gradually by the neutrino diffusion with the time 
scale of 10-100 seconds and evolves to a nearly “cold” NS 
with Yi, ~ 0 and S ~ 0, unless another collapse to a black 
hole does not take place [M1IM1I57]. 

Thermal properties of the hot NSs are intimately re¬ 
lated to the physics of high density matter at finite tem¬ 
perature. Indeed, the hot NSs with the hadron-quark mixed 
phase (Fig.j^a)) have been studied previously, e.g. [551IMI 
isniisi]. It is then interesting to explore properties of the 
hot NSs (such as the M — R relation and the profiles of 
the temperature, density and sound velocity inside the 
star etc) with hadron-quark crossover (Fig. [^b)). 

In Figj^ we show a schematic picture which compares 
the internal structure of the hot and cold NSs with 2Mq. 
Above the normal nuclear matter density po, we use the 
EOS interpolated between TNI2 (hadron) and NJL (quark). 
On the other hand, below po, we use the thermal EOS 
which consists of an ensemble of nuclei and interacting nu¬ 
cleons in nuclear statistical equilibrium given by Hempel 
and Schaffner-Bielich (HS EOS) [SI]. (Use of other EOSs 
below pq does not show quantitative difference as dis¬ 
cussed in [63] . ) Once the baryon density becomes smaller 


than the neutron drip density the temperature 

becomes smaller than O.I MeV. Then we switch to the 
standard BPS EOS |M| . 


7.1 Supernova Matter at finite T 

Major differences of the supernova matter from that of 
the cold matter discussed in previous sections are (i) the 
diffused Fermi surface due to the effect of T, (ii) the exis¬ 
tence of degenerate neutrinos, and (iii) the contributions 
from anti-particles. 

By neglecting the strangeness in hadronic matter and 
the muon which are irrelevant for stiff EOS, we consider a 
system composed of n, p, e~, e+, Ve and Pe in the hadronic 
EOS at finite T, and u, d,s, e~, e+, and Pg in the quark 
EOS at finite T. The charge neutrality, chemical equilib¬ 
rium and baryon and lepton number conservations are im¬ 
posed. In practice, we find Pe/T > 15 in the interior of the 
hot NSs, so that the effects of e+ and P^ (as well as neutri¬ 
nos in second and third generations) are negligibly small. 
The color superconductivity is switched off for simplicity. 

We use the Helmholtz free-energy per baryon F = 
F/N^ = E — TS as a basic quantity to interpolate the 
hadronic matter and the quark matter at finite T m 
This is a natural generalization of A at T = 0 in the pre¬ 
vious sections. E is a function of p, T and Yj, so that we 
have 

Fip,T-Yi)=FH{p,T-,Yi)w_ip,T) 

+Fq{p,T;Yi)w+{p,T). (14) 
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(a) (b) 

Fig. 10. (a) The temperature T of the isentropic matter as a function of the baryon density p for Yj = 0.3 with S' = 1 (red) and 
S = 2 (blue) in hot CRover EOS. The dashed line corresponds to the hot TNI2 EOS only with hadrons and leptons, (b) The 
isentropic pressure P of hot CRover EOS as a function of e for {Yi, S) =(0.3, 1), (0.3, 2) and (0.4, 1). The black line corresponds 
to the cold CRover EOS for cold neutron star matter. The crossover window is shown by the shaded area on the horizontal 
axis. These figures are adapted from mi. 




Fig. 11. (a) The neutron star mass M as a function of the total baryon number Nb in CRover EOS. Red and blue curves 

correspond to the hot and cold neutron stars, respectively, (b) The energy release AE = Mhot — Mcoid as a function of the cold 
neutron star mass Mcow. (Y, S) = (0.3,1) is adopted. These figures are adapted from [17] . 


Here Fp and Fq are the hadron-blepton free-energy per 
baryon and the quark-l-lepton free-energy per baryon, re¬ 
spectively. The typical temperature of hot NSs is about 
30MeV which is sufficiently smaller than the thermal dis¬ 
sociation temperature of hadrons of about 200 MeV. There¬ 
fore, we ignore the T-dependence of the weight functions, 
w±{p,T) w±{p). 

The entropy per baryon and the energy per baryon 
are obtained by using the thermodynamic relations, S = 


—dF/dT and E = F + TS. Under the assumption that 
w± is T-independent, Eq.(14) leads to 


S{p,T;Yi) = Su{p,T]Yi)w-{p) + S(^{p,T-,Yi)w+{p), 
E{p, T;Yi)= Fh [p, T; Yi)w_ (p) + Eq{ p, T; Yi)w+ (p). 


To obtain Fh in Eq.( 14), we solve the thermal Hartree- 
Eock equation for isothermal matter composed of n,p, e~, 
e"*", iZg and 9^ with the same effective nucleon interaction as 
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(C) (d) 

Fig. 12. (a) Mass-Radius (M-R) relationship for {Yi, S) = (0.3,1). Red: hot neutron stars with hot CRover EOS. Blue: cold 
neutron stars with cold CRover EOS. (b) The sound velocity squared Vg as a function of the distance from the center r of IAMq 
neutron star. Colors on each line are the same as in (a), (c) The density profiles of the hot neutron star with M = 1.4Mq and 
{Yi,S) = (0.3,1). Solid and dashed lines correspond to the EOS with crossover and without crossover, respectively, (d) The 
temperature profiles of the same neutron star as plotted in (c). These figures are adapted from |17| . 


the TNI2 and TNI2u EOS at T = 0 (details are shown in 
[S5]). To obtain Fq in Eq.( 14), we start with the Gibbs free 
energy calculated in the NJL model in 0 = 

f2quark(At, V,T] fii) + f7iepton(/i, V, T ; /T; jwith fi and fii be¬ 
ing the baryon chemical potential and the lepton chemical 
potential, respectively. Then we make the Legendre trans¬ 
formation from the Gibbs free energy Hq to the Helmholtz 
free-energy Fq [17]. 


In the following, we consider typical values of the lep¬ 
ton fraction Yj = 0.3, 0,4 and those of the entropy per 
baryon S = 1,2. The crossover window is fixed to be 
(p,F) = (SpoiPo); and = 0.5Gs is adopted. 


The thermodynamic quantities for isothermal matter 
with p, T and 1) as parameters can be converted into those 
for isentropic matter with p, S and Yi as parameters by 
using the relation T = T{p; S,Yi). In Fig, [^a), the tern- 
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Fig. 13. (a) The central temperature Tc and (b) the central density pcent as a function of the neutron star mass M of hot NS 
with (Vi, 5) = (0.3,1). The solid (dashed) lines correspond to the EOS with (without) crossover. These figures are adapted from 

m- 


perature T of the isentropic matter is shown as a function 
of the baryon density p for 1) = 0.3 with ,§ = 1 and S = 2. 

The isentropic pressure is obtained through the ther¬ 
modynamic relation 


P{p,T{p-S,Yi)-YuS)=p^ — 


(15) 


S,N,Yi 


The hot CRover EOS for supernova matter is shown in 
Fig. [^b) for {Yi,S) = (0.3,1), (0.3,2) and (0.4,1). For 


comparison, the cold CRover EOS for cold neutron star 
matter (T = 0 without neutrino degeneracy) is shown 
by the black solid lines. The hot CRover EOS and cold 
CRover EOS are similar except for the low density region 
where the hot EOS becomes stiffer. 


7.2 Properties of Hot NSs 


In Fig. [TT|(b), we show AE = Mjjot — M^oid in the unit 
of Mq as aOunction of the mass of cold NS, Mcom- The 
typical amount of energy released due to the contraction 
is about O.O4M0 for Mcom = 1.4M0. 

In Fig. 12 (a), we plot M — R relation for hot and cold 
NSs with and without crossover. The maximum mass of 
hot NSs is very similar with that of cold NSs. On the other 
hand, a big difference of the radius can be seen between 
hot and cold NSs. This comes from the stiffening of the 
hot EOS at densities below pg. The local sound velocity 
squared Vg{r) for isentropic matter can be calculated as 


'^S \P^ ^ 5 


dP{p,T{p);Yi,S)/dp 


Yi,s de{p,T{p);Yi,S)/dp 


(16) 


Yi,S 


Fig. 11 (a) shows the gravitational mass M as a function 
of the total baryon number Nb for hot (red line) and 
cold (blue line) neutron stars. The hot neutron stars have 
larger mass than the cold ones for given Nb- We note 
that the maximum value of Nb for hot NSs is smaller 
than that for cold NSs. This implies that hot NSs do not 
have a chance of delayed collapse into black holes after 
the cooling. This is in contrast to the previous works with 
the exotic phases such as the pion condensation [66] and 
the hadron-quark mixed phase HSIISO]; in those cases, the 
softening of the EOS due to exotic phases is tamed by the 
finite temperature effect, so that the maximum value of 
Nb for hot NS becomes larger than that of the cold NS 
and there arises a chance of the delayed collapse. 


with p(r) obtained by the TOV equation. In Figj^ (b), 
sound velocity squared is plotted as a function of tlie dis¬ 
tance from the center r for M = 1.4M0. The sound veloc¬ 
ity in cold NS is larger (smaller) at higher (lower) density 
than that of the hot NS [T7] . 

To see the thermal and lepton effects on the internal 
structure of the hot NSs, we plot the density profile p{r) 
and the temperature profile T{r) of hot NSs with canonical 
1.4M0 in Fig. 12 (c) and (d), respectively. Due to the 


stiffness of the CRover EOS, the central density becomes 
smaller and the profile becomes flatter as shown in Fig. 
[T^c). Moreover, as we have shown in Fig[l0j(a), T becomes 
smaller for given p by the crossover to quark degrees of 
freedom. Those are the reasons why the temperature of 
the internal core of the hot NS becomes smaller and flatter 
with crossover (Fig. [T^d)). 

In Fig. 13 (a) and (b), the central temperature Tcent 
and the central density Pcent of the hot neutron stars with 
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and without the crossover are plotted as a function of neu¬ 
tron star mass M. The effects of crossover on the internal 
structure of the NSs shown in Fig.[I^c,d) for M = IAMq 
can be seen for wide range of M. 


8 Summary and Concluding Remarks 

In this article, we have studied the bulk properties of 
cold and hot neutron stars on the basis of hadron-quark 
crossover picture. A new EOS, “CRover”, introduced in 
[T^fTHlfTT] describes the smooth transition from hadronic 
matter to quark matter in a phenomenological way. The 
hadron-quark crossover which makes the EOS stiffer by 
the effect of the quark matter is in contrast to the first- 
order hadron-quark transition leading to the softening of 
the EOS. 

At zero temperature, a crossover at around 3po leads 
to an EOS hard enough to sustain the 2Mq NSs. The radii 
of the NSs are located in a narrow region (12.5 ± 0.5) km 
which is insensitive to their masses. Due to the stiffening 
of the EOS induced by the crossover, the central density 
of the NSs is at most 4 pQ. Therefore, hyperon mixing, 
whose threshold density is about 4 po in the CRover EOS, 
takes place only for very massive NSs. This constitutes a 
solution of the long-standing hyperon puzzle. 

We have studied the effect of color superconductivity 
(CSC) on the bulk properties of NSs under the hadron- 
quark crossover. With the diquark coupling H/G^ = 1, we 
find that 2SC phase may appear inside the NSs while the 
onset of the CEL phase is too high to be realized even in 
massive NSs. As a result of a slight softening of the EOS 
due to the color superconductivity, the maximum mass of 
the NSs with CSC is reduced by about 0.2 Mq from that 
without CSC. 

To examine the properties of the hot NSs at birth at 
finite temperature, we considered the supernova matter 
with the CRover EOS generalized to isentropic environ¬ 
ment. We found that the hadron-quark crossover plays an 
important role to lower the central temperature of hot 
neutron stars in comparison to the case of hadronic EOS. 
This suppression of temperature comes from the presence 
of the quark degrees of freedom in the crossover region. 
Given baryon number, hot neutron stars have generally 
larger radius and larger gravitational mass caused by the 
high lepton fraction and the thermal effect. This suggests 
that, during the contraction from hot NS to cold NS, grav¬ 
itational energy is released and simultaneously the spin- 
up takes place. The released energy is shown to be about 
0.04 Mq and the spin-up rate is about 14 % (assuming the 
conservation of angular momentum) for M^oid = 1.4M0 
of evolved cold neutron stars. 

The hadron-quark crossover turns out to have inter¬ 
esting phenomenological implications to the key issues of 
the neutron stars, such as the massive neutron stars and 
hyperon puzzle, the universal radius of the neutron stars, 
temperature and density profiles inside the hot neutron 
stars, and so on. One of the most important and yet chal¬ 
lenging theoretical problems is to elucidate the QCD ba¬ 


sis of the phenomenological hadron-quark crossover intro¬ 
duced in this article. 
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